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Abstract We explore features of redshift distortion in Fourier analysis of N-body 
simulations. The phases of the Fourier modes of the dark matter density fluc- 
tuation are generally shifted by the peculiar motion along the line of sight, the 
fvq \ induced phase shift is stochastic and has probability distribution function (PDF) 

CN| ■ symmetric to the peak at zero shift while the exact shape depends on the wave 

^J \ vector, except on very large scales where phases are invariant by linear perturba- 

tion theory. Analysis of the phase shifts motivates our phenomenological models 
for the bispectrum in redshift space. Comparison with simulations shows that our 
toy models are very successful in modeling bispectrum of equilateral and isosceles 
triangles at large scales. In the second part we compare the monopole of the power 
spectrum and bispectrum in the radial and plane-parallel distortion to test the 
O \ plane-parallel approximation. We confirm the results of Scoccimarro (2000) that 

difference of power spectrum is at the level of 10%, in the reduced bispectrum such 
difference is as small as a few percents. However, on the plane perpendicular to 
the line of sight of k^ = 0, the difference in power spectrum between the radial 
and plane-parallel approximation can be more than ^ 10%, and even worse on 
very small scales. Such difference is prominent for bispectrum, especially for those 



5-H ' configurations of tilted triangles. The non-Gaussian signals under radial distortion 

on small scales are systematically biased downside than that in plane-parallel ap- 
proximation, while amplitudes of differences depend on the opening angle of the 
sample to the observer. The observation gives warning to the practice of using the 
power spectrum and bispectrum measured on the fc^ = plane as estimation of 
the real space statistics. 

Key words: cosmology: theory — large-scale structure of universe — methods: 
statistical 

1 INTRODUCTION 

Large redshift surveys of galaxies provide us with distances estimated from measured redshifts 
to complete our three-dimensional topography of the visible universe. However, due to peculiar 
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velocities of galaxies, the inferred redshift from observation is not an exact indication of the true 
distance, therefore our estimation of statistics in redshift space suffers from redshift distortion. 
On large scales peculiar velocities are dominated by bulk inflow, the redshift distortion boosts 
the clustering strength of galaxies, on small scales the virialised peculiar velocities suppress 
the clustering power and account for the phenomenon of "Finger of God" (e.g. Kaiser 1987; 
Hamilton 1998). 

Redshift distortion is entangled with bias and the gravitational nonlinearity. Nevertheless, 
fruitful results of redshift distortion have been achieved and applied in statistical analysis of 
LSS, particularly the power spectrum P{k) and the two point correlation function ^. On large 
scales where linear theory applies, the power spectrum in redshift space is 

P,{k) = {l + [3fi^fP{k), (1) 

where /3 = fl'^ /b with b the linear bias parameter, and /i is the cosine of the angle between k 
and the unit vector z of line-of-sight (LOS) (Kaiser 1987). Further observation and analysis of 
N-body simulation denotes that over wide scales including strong nonlinear regime there is an 
accurate empirical formula to model the effect of redshift distortion on power spectrum with 
introduction of a pairwise velocity dispersion parameter a^ (e.g. Park et al. 1994; Peacock & 
Dodds 1994; Cole, Fisher & Weinberg 1995; Hatton & Cole 1998; White 2001), 

''•« ^ ^'^' [ll'(W/2]' ^ ''' 

In comparison, understanding of redshift distortions for bispectrum or the three point cor- 
relation function is not satisfactory (Hivon et al. 1995; Verde et al. 1998; Scoccimarro ct al. 
1999), which, in part, is attributed to the non-perturbative nature of redshift distortion effect 
(Scoccimarro et al. 1999). The intrinsic difference between power spectrum and bispectrum is 
that power spectrum is a phase invariant quantity which contains only the information of the 
amplitude of the Fourier transform of the density contrast, bispectrum consists of contributions 
from both amplitudes and phases (Matsubara 2003). Research based on phases has revealed 
that phase information is very crucial in LSS distribution pattern recognition (Chiang & Coles 
2000; Chiang 2001; Chiang, Coles & Naselsky 2002), especially for the bispectrum and the three 
point correlation function (Watts & Coles 2003; Watts, Coles & Melott 2003; Matsubara 2003; 
Hikage, Matsubara & Suto 2004; Chiang 2004; Hikage et al. 2005). Therefore there is strong 
reason to believe that the behavior of phases under redshift distortion will provide us with 
precious information, which is the starting point of this paper. 

An alternative route to tackle redshift distortion is to recover real space quantities from 
measurements in redshift space. The anisotropic two point correlation function in redshift space 
^(s) is measured along two separations as ^(cr, tt), where tt = s • z/s and cr = |s x z\/s. The 
projected two point correlation function S(ct) is free of redshift distortion and simply related 
to the real space ^(r) by (Davis & Peebles 1983; Fisher et al. 1994) 

E{a) ^ r^ ^{a, n)dn = 2 /^°° ^^M^r^ ^ (3) 

J -00 J a (r^ - cr2)2 

or, we can deproject S to have the real space two point correlation function (e.g. Saunders, 
Rowan-Robinson & Lawrence 1992; Hawkins et al. 2003) 

e(r) = / ^-^^da . (4) 
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In Fourier space it is much simpler, as we can see from Eq. 1 and Eq. 2 that if we restrict our 
measurement on the plane of /c^ = only, what we obtain is the real space power spectrum 
(e.g. Park et al. 1994; Hamihon & Tegmark 2002). 

However, the approach relies on the assumption of the parallel approximation, i.e. the 
observer is at infinite distance so the unit vector z of the LOS is pointing to a fixed direction. 
In reality there is no such fixed direction of LOS, systematic biases are inevitably introduced in 
the treatment mentioned above. Particularly since Fourier transform mixes all scales together, 
even at high fc there is some contamination from scales at wide angles. Consequently not only 
are polyspectra measured on /c^ = plane not equivalent to those in real space, but also 
the monopoles of polyspectra with respect to LOS are not strictly following what theories of 
redshift distortion under plane approximation predict. Further complication also arises that 
Fourier modes are no longer independent and the covariance matrix is not diagonal (Zaroubi & 
Hoffman 1996, Scoccimarro 2000). 

Attention has been called in works analyzing galaxy redshift surveys which always subtend 
wide angles over sky and are at low redshift, to reduce the systematics in the adoption of 
parallel approximation. Numerical studies have shown that if one confines work with small 
angle, differences brought forward by parallel approximation in monopoles are negligible, in 
quadruples of power spectrum are usually under ~ 5% if the opening angle of the window in 
which Fourier transform is performed is less than 50° (e.g. Park et al. 1994; Cole, Fisher & 
Weinberg 1994). Scoccimarro (2000) also reported that such approximation was good enough 
for measuring the monopoles of the power spectrum and the bispectrum. Even for all sky 
surveys the difference is as small as 10%. Whatsoever, as the parallel approximation plays 
a fundamental role in Fourier analysis of galaxies' spatial distribution, it definitely deserves 
exclusive studies, which is very crucial in keeping confidence on our interpretation of what we 
estimated in redshift space. We will show in this paper that the real life is never simple. 

Our paper is organized as follows. In Section 2, we give general description of mapping 
from real space to redshift space, demonstrate the shift of phases of Fourier modes by redshift 
distortion and discuss possible construction of empirical model for bispectrum in redshift space. 
Section 3 focuses on the difference between the plane-parallel approximation and the radial 
distortion. Conclusion and discussion is in section 4. Throughout the paper we are using two sets 
of ACDM simulations by the Virgo consortium, the Virgo simulation in box of size 239.5ft.~^Mpc 
with 256^ particles and the VLS simulation in box of size 479/i~^Mpc with 512^ particles, 
cosmological parameters of the two simulations are 17^ = 0.3, f^A = 0.7, cs = 0.9 and F = 0.21 
(Jenkins et al. 1998). 

2 PHASE INFORMATION AND BISPECTRUM IN REDSHIFT SPACE 

2.1 Phase shift is a nonUnear effect 

For simplicity in this section we adopt the plane approximation so the unit vector of LOS z has 
fixed direction. The mapping from real space coordinates r to redshift space s is 

s = r-/3uz(r)z , (5) 

where u(r) = — v(r)/(7i/3), v(r) is the peculiar velocity, and H-ir) = [1 / a){da / dr) = Ha is the 
conformal Hubble parameter with FRW scale factor a and conformal time r (Scoccimarro et 
al. 1999; Scoccimarro 2004). 

The density contrast in redshift space, (5s (s), is obtained from the real space density fluc- 
tuation 5{r) by requiring conservation of the number of galaxies, 

{l + 5s)d^s = {l + 5)d'^r (6) 
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Fourier transform of the number density contrast in redshift space reads 

5,(k) ^ I ^e-^'^-S^is) = I ^e-^^-^e^^'^M^^ [<5(x) + /3V.u,(x)] , (7) 

in which and hereafter the subscript "s" refers to quantities in redshift space. The term in square 
brackets describes the squashing effect, i.e., the boost to the clustering ampHtude, whereas the 
exponential factor encodes the "Finger of God" effect, which erases power due to velocity dis- 
persion along the LOS. In plane-parallel approximation, the distortion comes from the peculiar 
velocity along the LOS Vz- Thus, to generate a redshift space sample, we only need s^ — rz — fiu^. 
It is clear from Eq. 7 we see that in the redshift space, not only the amplitude |(J(fc)|, but 
also the phase angle 9{k) is distorted which is casually ignored. The expression can be expanded 
perturbatively as demonstrated by Soccimarro et al. (1999) who calculated the expansion to 
second order to model the bispectrum in redshift space on large scales. We reproduce here the 
expansion 

5.(k) = V [ d^h...d^k4SD]n [(5(ki) +/3Ai?0(ki)] ^f^^^" ' ^0(k2)...^(/)(k„) , (8) 

~i^J (n-iy. k2 kn 

where [(Juln = <5i5(k — ki — ... — kn) is the Dirac function, velocity divergence is (f>{r) = V • u and 
/i„ = k„ • z/kn. In the linear theory only the first order n — 1 counts, we recover the well-known 
Kaiser formula (Kaiser 1987) 

6s{k)^{l + P^i^)S{k) , (9) 

which tells that at the first order phases are not changed by redshift distortion, phase shift is a 
phenomenon at high orders n > 2, i.e. nonlinear. We are not going to calculate quantitatively 
the phase shift in perturbation theory, rather, as a starting point in this direction, we aim at 
presenting the existence of the phase shift in numerical simulations and henceforth implication 
drawn from the information to understand redshift distortion on phase-related statistics, the 
bispectrum. 

2.2 Bispectrum and phase information 

By definition, the bispectrum is the ensemble average of the product of three Fourier modes 
with their wave vectors forming a triangle, 

B(ki,k2,k3--ki-k2) = (<5(ki)J(k2)<5*(k3)) ^ {Md2\\S3\e'^''+'''''^) , (10) 

where we wrote the Fourier mode S(k) — |(S|e*^. It is well known that the nonlinearity of 
gravitational evolution generates correlation between the phase 9 and the amplitude \S\ as well 
as the correlations between different modes, especially on small scales although in the initial 
fluctuation there are no such correlations. Whilst we are living with such universal correlations, 
on large scales where the nonlinearity is fairly weak, keeping the correlation between different 
modes to ensure non-zero bispectrum we can make an assumption to neglect the correlation 
between the amplitude and the phase of a single mode, an interesting conclusion emerges. With 
the assumption, the ensemble average of bispectrum can be decomposed into 

S(ki,k2,k3) cx ilSilMSsl) J e'''p{Q)dQ , (11) 

in which p{0) is the PDF of the phase sum Q = 9i + 62 — 03. Since the bispectrum is real 
and (5(k) is Fourier transform of a real function, p{0) shall be symmetric to O — and can 
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be expanded as p{0) ~ cq + '^^^ Cn cos(n6). Keeping only the first two terms, we recover the 
essence of tlie results of Matsubara (2003) 

p{e) (x{l + const. -cos e) . (12) 

In redshift space, modifications occur in both amplitudes and phases, 

<5,(k) = |,5,|e^'^= = a(k)|<5(k)|e'(^+^^) , (13) 

in which the amplification factor a(k, /i) is real and phase shift is A0(k) = ^^(k) — 0(k). Thus 
the bispectrum in redshift space can be written as 

B,(ki,k2,k3 = -ki -k2) = {a,a2as\Si\\S2\\Ss\e'^''+''-'^'^e'^^'^+^<'--^'^^^) . (14) 

According to our observation of the Virgo and the VLS simulations, a(k) = a{k, ii) is deter- 
ministic and can be well approximated by a = (1 + /3/i^)/[l + (kfiay)'^ /2] as in Eq. 2, hence 
the product of amplitude amplification factors in Eq. 14 can be placed outside the ensemble 
average. 

If we assume at a given k. Ad is fixed regardless of the original stochastics, AO := A^i + 
A02 — A03 will be a constant so that Bg ~ aia2a^Be'^'^^\ A8 has to be mr [n — 0, ±1, ±2...) 
as Eg is real, so 

Bg = aia2a^B . (15) 

It is hard to conjecture there is cosmic conspiracy to satisfy AO ~ m: for arbitrary (ki, k2, ka) 
triplet, the only possibility is for any mode A9 ~ mr, which is obviously impossible as implied 
by Eq. 7 & 8. The immediate conclusion is the phase shift resulted from redshift distortion has 
dependence on the phase in real space and therefore is stochastic. 

In N-body simulations we found that although phases no matter in real space or in redshift 
space are uniformly distributed, the distribution of phase shift p{A9) is not uniform as shown 
in Fig. (1). In general p{A6) depends on (fc,/i) and is peaked at center A6 = 0. With the 
decreasing of fi and fc, the dispersion of the distribution is becoming smaller and more modes 
concentrate on the vicinity of A9 = 0. This is easy to understand since those modes around the 
plane perpendicular to LOS suffers little from redshift distortion (Eq. 7) while on large scales 
the deviation to linear redshift distortion theory of Eq. 9 becomes small. We examined that in 
simulations phases of modes on the plane p — are unshifted. 

2.3 Toy models for bispectrum in redshift space 

In attempt to understand the redshift distortion of bispectrum, as a crude approximation in 
light of the spirit of deriving Eq. 12 we assume that AO is uncorrelated with |(5i,2.3| and 
^1,2,3, which is possible since AQ is a sum of three random variables. Eq. 14 then becomes 
Bg = B X aia2as J e"^^^ p{AQ)dA<d , and p{AQ) ex (1 + const. ■ cos AG) if high order terms are 
ignored. With a further assumption that A0i,2,3 is independent of each other, we have 

Bs(ki,k2,k3) =B(ki,k2,k3) x aia2a3 

X f e''^^'p{A0i)dAei f e'^^^p{A92)dA02 j e'^'^^p(Ae:i)dAes. , ^ ' 

from which trivially we again have features of p{A6) that it is symmetric to A6 = seen in 
Fig. 1 and ex (1 + const. ■ cos A9). Note that here a and A9 are functions of two variables (fc, n). 
Eq. 16 is quite similar to the phenomenological model of Scoccimarro et al. (1999), and it serves 
as a theoretical support to their ansatz. 
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Fig. 1 Probability distribution functions of phase shift A6' resulted from redshift 
distortion in N-body simulations. The left panel plots the p{A9) on different k scales at 
fixed /i = k-z/fc, and in the right panel we show p{A9) at variant ^ at k = 0.4/iMpc~-'^. 



In practice often it is simpler to measure of the bispectrum averaged over all possible 
combinations of wave vectors (ki,k2,k3 = — ki — k2) with fixed fci — |ki|, ^2 — |k2| and the 
angle ip = ki-k2/(fcifc2). In real space it is the isotropic bispectrum B{ki, fc2, /ca) = -B(ki, k2, ks), 
in redshift space due to the anisotropy of redshift distortion, the angular averaged with respect 
to LOS is the monopole of true bispectrum _Bs(fci, A:2,'0iM1jM2jM3)- In that we denote the 
azimuthal angle of k2 about ki with Lp, and 



/xi = /i = ki • z/ki , 

M2 == M cos -0 — yl — /i^ sin ip cos ip , 



M3 



ki fc2 



(17) 



the monopole of bispectrum in redshift space is 



B(")(fci,fc2»-^ r 



dn / dipBs{ki,k2,xp,^J.,(p) . 
1 Jo 



(18) 



On large scales where the nonlinearity is weak and phases shifts are close to zero (Eq. 9 and 
Fig. 1), we can approximate the modulation by peculiar velocity with Eq. 15, 
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(19) 



We measured angular averaged power spectrum and bispectrum both in real space and 
redshift space for the Virgo simulation and eight subsets of the VLS simulation which are 
generated by dividing the VLS simulation into independent cubes of half of the original box 
size. (Tt, = 4.87 is fitted from power spectra with Eq. 2, and then is inserted into Eq. 19 as our 
toy model for i?s . In Fig. 2 it is interesting to observe the remarkable agreement of our toy 
model based on phase information with simulations for equilateral triangles for k < 0.2/iMpc~^ 
beyond which tree-level perturbation theory and the second order Lagrangian theory also breaks 
(Scoccimarro 2000). As k grows, phase shift is no longer negligible and our toy model fails. 
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Fig. 2 Bs /B of equilateral triangles. Dotted line is the prediction of Eq. 19, and 
dash line is the empirical model of Eq. 20. Points are averaged over measurement of 
the Virgo simulation and the eight subsets of the VLS simulation. Error bars are the 
Icr scatter. 



As a next level advanced than the simple prescription of Eq. 19, a better approximation 
may be achieved through integration of Eq. 16 instead of the simple formula of Eq. 19 if we 
know p(A0) which carries important information about the redshift distortion and numerically 
contains more dependence on (fc, /i) (Fig. 1). Because we arc not equipped with a whole battery 
of many simulations in very large box, so far in this work we are not developing functional 
expression to fit p(A0)'s dependence on (fc,/i) and a^. Nevertheless, we can parameterize the 
toy model in a simple way to improve its performance just as in Scoccimarro et al. (1999), e.g. 
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(20) 



with two fitting parameters Ai_2- An eye-ball examination with Ai = 1/1.15 and A2 — 1.75 
does improve agreement up to scale k ~ 0.5/iMpc~^ for equilateral triangles (Fig. 2). We then 
compare the modified toy model for other configurations with simulations in Fig. 3. We see 
that for isosceles with ip >^ 90° the agreement of our modified toy model with simulations is 
very impressive too, although for those very tilted configuration with ^2/^1 > 2 differences are 
huge. We are not intending to give accurate fit, we are just demonstrating the possibility of 
constructing a good approximation model with the simple prescription of Eq. 19, and, in the 
mean time showing the discrepancy due to our rough handling with phase shift. 



3 GOODNESS OF PLANE-PARALLEL APPROXIMATION TO 
RADIAL DISTORTION 

So far we are working with plane-parallel approximation. As we have addressed in section 1, 
such approximation is appropriate when the angle subtended by the sample is small. When 
dealing with large samples, especially all sky surveys, radial distortion should be taken into 
account. To investigate the goodness of the plane-parallel approximation to the exact radial 
distortion, we are to compare a set of power spectra and bispectra in redshift space in the 
parallel approximation and the radial distortion. 
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Fig. 3 Bs /B{ki,k2,tp) of different configurations. The dotted line is our modified 
toy model Eq. 20, see Fig. 2 for other keys. 
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3.1 Monopoles of Power Spectrum and Bispectrum 

It has been shown by Scoccimarro (2000) that the monopoles of power spectrum and bispectrum 
of an all sky survey, in the radial and the plane-parallel distortions are consistent with each 
other within error bars. For the power spectrum the difference is within ^ 10% and for the 
bispectrum is within a few percents. 

To generate samples in rcdshift space in radial distortion, we strictly follow Eq. 5 with 
z defined by the unit vector pointing from the observer to the point of sample, i.e. z is 
simply the unit position vector f. Cartesian coordinates ^^=1^2, 3 of any point in our sam- 
ples are within [0, L]/i~^Mpc, L — 239.5. In this subsection, we placed the observer at 
(0.5\/2i, 0.5\/2L, — 100)/i^^Mpc, the largest angle of the sample opening to the observer is 
118.88°. Then points in samples are transformed into the coordinate system in which the ob- 
server is sitting at the origin point. For every data point we have redshift from Hubble fiow zuub 
and peculiar velocity induced Zpec, the comoving distance in rcdshift space s is then calculated 
from the final redshift z — {1 + ZHub){^ + Zpec) — 1 while its orientation docs not change. 

In Fig. 4, we plotted the ratio of the monopoles of power spectra in radial distortion and 
the parallel approximation, P} / Pp where the subscript "r" means radial distortion, "p" 
refers to parallel approximation and all polyspectra are in redshift space if not specified here- 
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Fig. 4 Ratio of the monopolcs of power spectra in radial distortion and parallel 
approximation. The sample is at the distance of 100/i^^Mpc. 



after. Monopole of power spectrum is defined as the angular average of the anisotropic -P(k). 
Comparison of the monopoles of reduced bispectrum is illustrated in Fig. 5. Note that the 
monopole of reduce bispectrum is given by 



Q'-"Hki,k2,ib) = 



B(o)(fci,fc2,V') 



P(0)(fci)P(")(A:2) + P^°Kk2)P^°Hh) + P(o)(fc3)P(°)(fci) 



(21) 



We can see from Fig. 4 that the differences between monopoles of power spectrum in the 
radial and the plane-parallel distortion are within 10%, while the differences for bispectra are 
of a few percents (Fig. 5). These results confirm the assertion of Zaroubi & Hoffman (1996) and 
Scoccimarro (2000). As shown in Fig. 4, the power in parallel approximation on large scales 
k < O.l/iMpc"^ is smaller than that under radial distortion. On the other hand, on small scales 
k > 0.2ft,Mpc~^, the plane-parallel distorted power spectrum is larger. This is because in the 
plane-parallel approximation there is no distortion effect along directions perpendicular to the 
fixed LOS, the actual redshift distortion is radial and the LOS is not fixed. Distortion is more 
severe to structure clumps so that there are more effects of boost at large scales and suppression 
at small scales on the clustering strength. 



3.2 Po^ver Spectrum and Bispectrum on the Plane of kz 







We have seen that there is systematic bias caused by the parallel approximation to the radial 
distortion in monopoles of power spectrum and bispectrum though very weak. In this subsection 
we set about to explore the behavior of the polyspectra on the plane of fcz =0 which are 
regarded as estimation of the real space quantities in the plane-parallel approximation. It is 
easy to understand that on large scales where the opening angle to the observer is wide so 
the lines of sight can not be approximated by parallel lines, such deviation is not separable 
in Fourier transform which is a global operation. Only when the opening angle of the sample 
window is small we are safe. In order to check the influence of the opening angle of sample 
to the method of using polyspectra oi kz = as estimation in real space, we created different 
samples in redshift space from simulations by placing the observer at distances of 100, 200 and 
1000/i~^Mpc respectively. 
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Fig. 5 Ql°^/Q'p*\ki,k2,^p) for triangles of h ^ 0.2, 0.3, 0.4, 0.5, O.e/iMpc^i and 
^2/^1 ~ 1, 2, 3. The observer is at the distance of 100/i^^Mpc. 



Fig. 6 shows that when k is larger than 0.2/iMpc^ , under the plane-parallel approxima- 
tion, Pp can grows from ^ 10% to more than 50% larger than P^; when k is smaller than 
0.2ft,Mpc~^, Pp could still be > 10% smaller than Pr- Meanwhile, as we expected, with incre- 
ment of the distance of the sample to observer, i.e. the opening angle getting smaller, such 
difference systematically decrease. 

However, even when the distance is 1/i^^Gpc which corresponds to an opening angle as 
small as ~ 20°, at fc = 3/iMpc^^ the systematical bias persists around 50% if one is interested 
in the highly nonlinear regime, meanwhile the bias is also > 10% for k < 0.03ft.Mpc~^ which 
prevents us from having measurements in precision. It is a bad news for most redshift surveys, 
such as the 2dFGRS redshift survey and the SDSS of which subtended angles are much larger 
than 20°. Furthermore, as seen in Fig. 7 and Fig. 8, there are similar systematic biases in 
reduced bispectrum measured at kz = to the real space bispectrum if we adopt the plane- 
parallel approximation, the biases not only attenuate the amplitudes but also the configuration 
dependence. 
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Fig. 6 Ratios of power spectra for radial and parallel distortion on the plane of k^ — 0. 
The solid line, dash line and dot line are for samples at distances of 100/iMpc~^, 
200/iMpc^^ and lOOO/iMpc^^ to the observer respectively. 



Therefore, it is hopeless to accurately estimate the real space polyspectra on large scales 
with the convenient method of using what we measured on the plane A:^ = in redshift space 
if no other special skills are applied to compensate the systematical bias. One of the possible 
application of the method is to estimate small scale polyspectra in real space by splitting the 
wide angle sample into many patches of very small angles in compromise with the scales of 
interests, and later collecting measurements in these small patches for analysis in precision, 
which, of course, needs careful numerical experiments to decide the angular widths of those 
patches (e.g. Park et al. 1994). 



4 SUMMARY AND DISCUSSION 

In the exploration of the characteristics of redshift distortion in Fourier space, we begin with 
research on phases of Fourier modes which are rarely noticed in previous works. On very large 
scales where linear theory applies, phases are not changed by redshift distortion, we conclude 
that phase shift Ad is a nonlinear effect. Then we find that the phase shift by redshift distortion 
is not deterministic, it has a distribution function p{A9) symmetric to its peak at A9 = 0, 
and qualitatively proportional to (1 + const. ■ cosA0). The distribution function has strong 
dependence on the wave length k and the orientation fj, of the wave vector to the LOS. More 
modes are populated around A9 — for small k and small /i. 

Limited by the sizes of the available Virgo and VLS simulations and computing resources, 
we can not obtain reliable p(A6) on scales of A; < 0.1/iMpc~^ to build empirical model for 
the distribution function. Although in principle it is possible to calculate the phase shift in 
perturbation theory, it is unclear how to generate the PDF of the phase shift in perturbation 
approach, which we leave to future work. Instead we make assumptions on the phase shift at 
different levels in effort to depict the redshift distortion on bispectrum. We discovered that 
our very crude consideration does provide us with a toy model (Eq. 19) for the monopole of 
bispectrum in redshift space which is in good agreement with numerical simulations on scales 
k < 0.2ft,Mpc"^, and probably can serve as a basis for extended phenomenological models 
such as Eq. 20 and the model in Scoccimarro et al. (1999). If in future an approximation for 
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Fig. 7 Ratios of reduced bispectra for radial and parallel distortion on the plane of 
kz = 0. Observer is located at the distance of lOOft-^^Mpc. 
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p{A9) analogous to the modeling of amplitude change like Eq. (2) can be set with help of large 
simulations, it is promising that integration of Eq. 16 will offer us with greatly improved model 
for the monopole of bispectrum in redshift space. 

On non-linear scales, the plane-wave components of density fluctuation do not evolve inde- 
pendently, interaction among difference modes leads to generation of coupled phases. Therefore, 
significant phase correlation between different modes is expected. As the Fourier phases of a 
random Gaussian field are randomly distributed, phase correlation, if any, should characterizes 
the non-Gaussianity of density fields (Matsubara 2003). Such correlation would have important 
effects on the bispectrum, which is the lowest order statistics of non-Gaussianity. The impor- 
tance of connection between phase correlation and non-Gaussianity is doubtless. By using a 
perturbative quadratic model. Watts & Coles (2003) have shown for the first time how phase 
association give rise to non-vanishing bispectrum and three point function in nonlinear pro- 
cesses. The arrangement of phases of modes determines weather such non-Gaussian descriptors 
are zero or not, but the magnitude of non-Gaussianity is determined by the magnitude of Fourier 
modes (Watts & Coles 2003). A relationship between phase correlations and the hierarchy of 
polyspectra in the Fourier space has been addressed by Matsubara (2003). Also, a numerical 
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Fig. 8 Ratios of reduced bispectra for radial and parallel distortion on the plane of 
kz = 0. Observer is located at the distance of 1000/i^^Mpc. 
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test on the bispectrum has been given. However, many details of this issue such as how phase 
correlation functions with redshift distortion and bias, remain unknown. 

The connection between polyspectra and phase information may render us much more 
useful information, while our knowledge of it is till limited. It is worthwhile for us to tackle the 
problem in more details by taking into account of both the phase shift of individual Fourier 
mode which is presented in this paper, and the correlation between phase shifts of different 
modes which importance is of priority in understanding the redshift distortion on bispectrum 
in strong nonlinear regime. As a first step, the simplest consideration on phase shift already has 
guided us in producing a good model, future inclusion of the PDF and the correlation between 
phase shift will shed light on the redshift distortion in Fourier space. 

In the second part of our report, we explored the goodness of the plane-parallel approxima- 
tion to the actual radial redshift distortion in Fourier analysis. The first examination is exerted 
to the monopoles of power spectrum and bispectrum. Our measurements with the Virgo and the 
VLS simulation are consistent with the results of Scoccimarro et al. (2000) that such deviation 
due to adoption of parallel approximation is ^ 10% for the monopole of power spectrum and a 
few percents for bispectrum. Note that the curve in Fig. 4 shows a systematic trend toward high 
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k. Then, we moved to check the rehabiUty of the approach which uses polyspectra measured 
on kz = plane as estimation in real space, which is true in plane-parallel approximation. We 
find that plane-parallel approximation to radial distortion brings in serious statistics systematic 
biases in this method. On large scales the "real space" power spectrum from measurements on 
fc^ = plane can be overestimated for more than 10%, and it is hugely underestimated more 
than 50% on small scales k >^ 3ft.Mpc~^ (Fig. 6). We also find that besides the amplitude of 
reduced bispectrum changes, the configuration dependence suffers from the badness of paral- 
lel approximation too (Fig. 8 & 9). According to our simulation, such deviation is persistent 
even the opening angle of the sample to the observer is as small as 20°. These discoveries pose 
questions on those claims based on the "real space" polyspectra obtained in this way. 

So as long as we are working with monopoles of power spectrum and bispectrum, it is secure 
to model redshift distortion with plane-parallel approximation. While if one is going to estimate 
real space power spectrum and bispectrum with measurements on plane kz = 0, one has to be 
careful to consider the caveats of the parallel approximation. Actually a more natural treatment 
is to decomposed the density contrast in redshift space with spherical harmonics and spherical 
Bessel functions (Heavens & Taylor 1995; Ballinger, Heavens & Taylor 1995; Fisher et al. 1995; 
Szalay, Matsubara & Landy 1998; Percival et al. 2004; Percival 2004; Szapudi 2004). 
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